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ABSTRACT 

We present a new model for the gas dynamics in the galactic disk inside the Sun's orbit. 
Quasi-equilibrium flow solutions are determined in the gravitational potential of the 
deprojected COBE NIR bar and disk, complemented by a central cusp and, in some 
models, an outer halo. These models generically lead to four-armed spiral structure 
between corotation of the bar and the solar circle; their large-scale morphology is not 
sensitive to the precise value of the bar's pattern speed, to the orientation of the bar 
with respect to the observer, and to whether or not the spiral arms carry mass. 

Our best model provides a coherent interpretation of many observed gas dynamical 
features. Its four-armed spiral structure outside corotation reproduces quantitatively 
the directions to the five main spiral arm tangents at |Z| < 60° observed in a variety of 
tracers. The 3-kpc-arm is identified with one of the model arms emanating from the 
ends of the bar, extending into the corotation region. The model features an inner gas 
disk with a cusped orbit shock transition to an X2 orbit disk of radius R ~ 150 pc. 

The bar's corotation radius is fairly well-constrained at R c ~ 3.5 ± 0.5 kpc. The 
best value for the orientation angle of the bar is probably 20 — 25°, but the uncertainty 
is large since no detailed quantitative fit to all features in the observed (I, v) diagrams 
is yet possible. The Galactic terminal velocity curve from HI and CO observations 
out to I ~ ±45° (~ 5 kpc) is approximately described by a maximal disk model with 
constant mass-to-light ratio for the NIR bulge and disk. 

Key words: Galaxy: structure, Galaxy: kinematics and dynamics, Galaxy: centre, 
Galaxies: spiral, Interstellar medium: kinematics and dynamics, Hydrodynamics. 



1 INTRODUCTION 

Although the Milky Way is in many ways the best stud- 
ied example of a disk galaxy, it has proven exceedingly dif- 
ficult to reliably determine its large-scale properties, such 
as the overall morphology, the structural parameters of the 
main components, the spiral arm pattern, and the shape of 
the Galactic rotation curve. A large part of this difficulty 
is due to distance ambiguities and to the unfortunate lo- 
cation of the solar system within the Galactic dust layer, 
which obscures the stellar components of the Galaxy in the 
optical wavebands. With the advent of comprehensive near- 
infrared observations by the COBE/DIRBE satellite and 
other ground- and space-based experiments, this situation 
has improved dramatically. These data offer a new route to 
mapping out the Galaxy's stellar components, to connect- 
ing their gravitational potential with the available gas and 
stellar kinematic observations, and thereby to understand- 
ing the large-scale structure and dynamics of the Milky Way 
Galaxy. 

From radio and mm-observations it has long been 
known that the atomic and molecular gas in the inner 



Galaxy does not move quietly on circular orbits: "forbidden" 
and non-circular motions in excess of 100 km s -1 are seen 
in longitude-velocity (I, w)-diagrams (e.g., Burton & Liszt 
1978, Dame et al. 1987, Bally et al. 1987). Some of the more 
prominent features indicating non-circular motions are the 
3 kpc-arm, the molecular parallelogram ( "expanding molec- 
ular ring" ) , and the unusually high central peak in the termi- 
nal velocity curve at I ~ ±2° . Many papers in the past have 
suggested that these forbidden velocities are best explained 
if one assumes that the gas moves on elliptical orbits in a 
barred gravitational potential (Peters 1975, Cohen & Few 
1976, Liszt & Burton 1980, Gerhard & Vietri 1986, Mulder 
& Liem 1986, Binney et al. 1991, Wada et al. 1994). 

In the past few years, independent evidence for a bar 
in the inner Galaxy has been mounting from NIR photome- 
try (Blitz & Spergel 1991, Weiland et al. 1994, Dwek et al. 
1995), from IRAS and clump giant source counts (Nakada 
et al. 1991, Whitelock & Catchpole 1992, Nikolaev & Wein- 
berg 1997, Stanek et al. 1997), from the measured large mi- 
crolensing optical depth towards the bulge (Paczynski et al. 
1994, Zhao, Rich & Spergel 1996) and possibly also from 
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stellar kinematics (Zhao, Spergel & Rich 1994). See Gerhard 
(1996) and Kuijken (1996) for recent reviews. 

The currently best models for the distribution of old 
stars in the inner Galaxy are based on the NIR data from 
the DIRBE experiment on COBE. Because extinction is im- 
portant towards the Galactic nuclear bulge even at 2 fim, 
the DIRBE data must first be corrected (or 'cleaned') for 
the effects of extinction. This has been done by Spergel, 
Malhotra & Blitz (1996), using a fully three-dimensional 
model of the dust distribution (see also Freudenreich 1998) . 
Binney, Gerhard & Spergel (1997, hereafter BGS) used a 
Richardson-Lucy algorithm to fit a non-parametric model 
of j(r) to the cleaned data of Spergel et al. under the as- 
sumption of eight-fold (triaxial) symmetry with respect to 
three mutually orthogonal planes. When the orientation of 
the symmetry planes is fixed, the recovered emissivity j(r) 
appears to be essentialy unique (see also Binney & Gerhard 
1996, Bissantz et al. 1997), but physical models matching 
the DIRBE data can be found for a range of bar orienta- 
tion angles, 15° < ipbect ?S 35° (BGS). (^bar measures the 
angle in the Galactic plane between the bar's major axis at 
/ > and the Sun-centre line. For the favoured iphm = 20°, 
the deprojected luminosity distribution shows an elongated 
bulge with axis ratios 10:6:4 and semi-major axis ~ 2 kpc, 
surrounded by an elliptical disk that extends to ~ 3.5 kpc 
on the major axis and ~ 2 kpc on the minor axis. 

Outside the bar, the NIR luminosity distribution shows 
a maximum in the emissivity ~ 3 kpc down the minor axis, 
which corresponds to the ring-like structure discussed by 
Kent, Dame & Fazio (1991), and which might well be due to 
incorrectly deprojected strong spiral arms. From the study 
of HII-Regions, molecular clouds and the Galactic magnetic 
field it appears that the Milky Way may have four main 
spiral arms (Georgelin & Georgelin 1976, Caswell & Heynes 
1987; Sanders, Scoville & Solomon 1985, Grabelsky et al. 
1988; Vallee 1995). These are located outside the 3kpc-arm 
in the {l,v) diagram and are probably related to the so- 
called molecular ring (Dame 1993), although it is unclear 
precisely how. One problem with this is that the distances 
to the tracers used to map out the spiral arms are usually 
computed on the basis of a circular gas flow model. The er- 
rors due to this assumption are not likely to be large, but 
cannot be reliably assessed until more realistic gas flow mod- 
els including non-circular motions are available. 

The goal of this paper is to construct gas-dynamical 
models for the inner Milky Way that connect the Galactic 
bar/bulge and disk as observed in the COBE NIR luminos- 
ity distribution with the kinematic observations of HI and 
molecular gas in the (l,v) diagram. In this way we hope to 
constrain parameters like the orientation, mass, and pattern 
speed of the bar, and to reach a qualitative understanding 
of the Galactic spiral arms and other main features in the 
observed (l,v) diagrams. 

In barred potentials, gas far from resonances settles on 
periodic orbits such as those of the x\- and X2-orbit families, 
and some important aspects of the gas flow can be under- 
stood by considering the closed periodic orbits (e.g., Binney 
et al. 1991). However, near transitions of the gas between or- 
bit families, along the leading edges of the bar, and in spiral 
arms, shocks form in the gas flow which can only be studied 
by gas dynamical simulations (e.g., Roberts, van Albada & 



Huntley 1979, Athanassoula 1992b, Englmaier & Gerhard 
1997). 

In this paper, we use the Smoothed Particle Hydrody- 
namics (SPH) method to study the gas flow in the grav- 
itational potential of the non-parametrically deprojected 
COBE/DIRBE light distribution of Binney, Gerhard & 
Spergel (1997), assuming a constant mass to NIR luminosity 
ratio. The gas settles to an approximately quasi-stationary 
flow, and the resulting model (I, v) diagrams enable us to un- 
derstand many aspects of the observations of HI and molec- 
ular gas. 

The paper is organized as follows. In Section 2 we give a 
brief review of the main observational constraints, followed 
in Section 3 by a description of our models for the mass dis- 
tribution and gravitational potential, and for the treatment 
of the gas. Section 4 describes the results from a sequence 
of gas dynamical models designed to constrain the most im- 
portant parameters, and compares these models with ob- 
servations. Our results and conclusions are summarized in 
Section 5. 



2 SUMMARY OF OBSERVATIONAL 
CONSTRAINTS 

2.1 Surveys and (l,v) diagrams 

The kinematics and distribution of gas in the Milky Way 
disk have been mapped in 21-cm neutral hydrogen emission 
and in mm-line emission from molecular species, most no- 
tably 12 CO. HI surveys include Burton (1985), Kerr et al. 

(1986) , Stark et al. (1992) and Hartmann & Burton (1997), 
surveys in 12 CO include Sanders et al. (1986), Dame et al. 

(1987) and Stark et al. (1988), and the most comprehensive 
survey in 13 CO is that of Bally et al. (1987). The large scale 
morphology of the gas based on these data is discussed in 
the review by Burton (1992). 

These surveys show a complicated distribution of gas 
in longitude-latitude-radial velocity (l,b,v r ) space or, inte- 
grating over some range in latitude, in the so-called (l,v) 
diagram. Fig. [j] shows an (l,v) diagram observed in 12 CO 
by Dame et al. (1998, see also Dame et al. 1987) and trans- 
formed by them to the LSR frame by subtracting \vq\ = 
20kms" 1 towards (l,b) = (56.2°, 22.8°) for the motion of 
the Sun with respect to the LSR (This corresponds to in- 
ward radial and forward components of the solar motion 
of uq = — 10.3kms _1 and vq — 15.3 kms -1 .) The general 
morphology of the (l,v) diagram in 12 CO is broadly similar 
to that obtained from HI 21-cm emission (see, e.g., Burton 
1992). 

In general, no additional distance information for the 
gas is available; thus - unlike in the case of external galaxies 
- the spatial distribution of galactic gas cannot be directly 
inferred. Converting line-of-sight velocities to distances, on 
the other hand, requires a model of the gas flow. 



2.2 Interpretation of (l,v) diagrams 

For the comparison of observed and model {l,v) diagrams, 
it is useful to first consider an axisymmetric disk with gas 
in circular rotation (e.g., Mihalas & Binney 1981). In this 
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Figure 1. (I, v) diagram for 12 CO from unpublished data by Dame et al. (f998). This figure contains all emission integrated over latitudes 
between b = —2° and 6 = 2°. The grey scale is adjusted such as to emphasize spiral arm structures. 



model, an observer on a circular orbit will find the following 
results: 

(i) Gas on the same circular orbit as the observer will 
have zero relative radial velocity. 

(ii) For clouds on a different circular orbit with velocity 
v(R), the measured radial velocity is 

v r = (w — u>o)Ro sinZ, (1) 

where u>(R) = v(R)/R is the angular rotation rate, Rq is 
the galactocentric radius of the observer, and the index 
to a function or variable denotes its value at Rq- From this 
equation we see that, as long as u)(R) decreases outwards, 
the radial velocities have the same sign as sin Z for gas inside 
the observer, and the opposite sign for gas on circular orbits 
outside the observer. 

(iii) For clouds inside the observers orbit (—90 < I < 
90°), the maximal radial velocity along a given line-of-sight 
Z is v r = sign(l)v(R) — VbsinZ. This so-called terminal ve- 
locity is reached at the tangent point to the circular orbit 
with R = Ro sinZ. For Z > (Z < 0) v r increases (decreases) 
from zero for clouds near the observer to the terminal ve- 
locity at a distance corresponding to the circular orbit with 
R — Rosinl; it then decreases (increases) again with dis- 



tance from the observer and changes sign when crossing the 
observer's radius at the far side of the Galaxy. The termi- 
nal velocities define the upper (lower) envelope in the (l,v) 
diagram for < I < 90° (-90° < I < 0). 

(iv) A circular orbit of given radius follows a sinusoidal 
path in the (l,v) diagram[cf. equation (1)], within a longi- 
tude range bounded by I — ± arcsin(i?/i?o). In the inner 
Galaxy (—45° < Z < 45°) circular orbits thus approximately 
trace out straight lines through the orgin in the (Z, v) dia- 
gram. 

(v) The edge of the Galaxy (the outermost circular 
orbit) results in a sine shaped envelope in the (Z, v) dia- 
gram with negative radial velocities for positive longitudes 
(0° < Z < 180°) and vice versa. The sign of v r on this enve- 
lope is different from that on the terminal velocity envelope. 

(vi) For a circular orbit model, one can derive the ro- 
tation curve of the inner Galaxy from the observed termi- 
nal velocities and eq. ([[]). This requires knowledge of the 
galactocentric radius Ro and rotation velocity Vo of the lo- 
cal standard of rest (LSR) as well as the motion of the Sun 
with respect to the LSR. 

Fig. [j] and similar HI 21 cm (l,v) diagrams (see, e.g., 
Burton 1992) show gas between 90° and —90° longitude 
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whose radial velocities have the wrong sign for it to be gas 
on circular orbits inside the solar radius. Yet this gas is ev- 
idently connected to other gas in the inner Galaxy and is 
not associated with gas from outside the solar radius. These 
so-called forbidden velocities, up to v r ~ 100 kms^ 1 , are a 
direct signature of non-circular orbits in the inner Galaxy, 
and they have been the basis of previous interpretations of 
the inner Galaxy gas flow in terms of a rotating bar poten- 
tial (e.g., Peters 1975, Mulder & Liem 1986, Binney et al. 
1991). 

One of the most prominent such features is the so- 
called 3 kpc-arm, visible in Fig. |l| as the dense ridge of emis- 
sion extending from (I ~ 10°, v — 0) through (I — 0,v ~ 
-50kms _1 ) to (I ~ -22°, v ~ -120 kms -1 ). 

At longitudes \l\ > 25°, the circular orbit model is a 
reasonable description of the observed gas kinematics. Most 
of the emission in 12 CO in fact comes from a gas annulus be- 
tween about 4 and 7 kpc galactocentric radius, the so-called 
molecular ring (e.g., Dame 1993), which probably consists 
of two pairs of tightly wound spiral arms (see Section 4) . A 
detailed interpretation of spiral arms in the (I, v) diagram 
requires a full gas dynamical model, as high intensities in 
Fig. |l] can be due to either high intrinsic gas densities or 
due to velocity crowding. 



2.3 Terminal velocities 

The top left and bottom right envelopes on the (l,v) dia- 
gram in Fig. |l| mark the terminal velocities. The terminal 
velocity curves will be used below for comparing with dif- 
ferent models and for cali brat ing the models' mass-to-light 

below, the terminal velocities 



300 



100 



ratios. As will be seen in §4.4 
in 12 CO and HI 21cm and between different surveys agree 
to a precision of ~ lOkms" 1 in most places, but there are 
some regions with larger discrepancies. 

Of particular interest is the strong peak in the terminal 
velocity curve with vt — 260 km s -1 at I ~ 2°. Outwards 
from there the drop in vt is very rapid; for a constant mass- 
to-infrared luminosity ratio it would be nearly Keplerian 
and would be hard to reproduce in an axisymmetric bulge 
model (Kent 1992). Instead, the rapid drop is probably con- 
nected with a change of orbit shape over this region (Ger- 
hard & Vietri 1986). In the model of Binney et al. (1991), 
the peak is associated with the cusped orbit in the rotat- 
ing barred potential, and the subsequent drop of v t with I 
reflects the shapes of the adjacent xi-orbits. 



2.4 Spiral arms 

From distant galaxies we know that spiral arms are traced 
by molecular gas emission. Indeed one can identify some of 
the dense emission ridges in Fig. [j] with Galactic spiral arms; 
where these meet the terminal velocity curve, they can be 
recognized as 'bumps' where dvt/dl ~ 0. In addition, spiral 
arms are clearly visible in the distribution of various tracers, 
such as HII regions. 

Fig. |^ shows an (I, v) diagram of several classes of ob- 
jects which are useful as discrete tracers of dense gas in spiral 
arms. On each side, we can identify two spiral arm tangents 
at around ± ~ 30 and ± ~ 50°. On the northern side, the 
~ 30° component splits up into two components at ~ 30° 




-100 - 



-300 



Figure 2. (I, v) diagram for HII regions from Gcorgelin & 
Georgelin (1976), Dowries et al. (1980) and Caswell & Haynes 
(1987) (open circles); and for massive molecular clouds from 
Dame et al. (1986) and Bronfman (1989) (filled circles). The ar- 
rows point to the positions of the densest clusters of warm CO 
clouds along the terminal curve, which are presumably located in 
spiral arm shocks (Solomon, Sanders & Rivolo 1985). For the sake 
of clarity, clouds with less than 10 5 ' 5 Mq have been omitted from 
Bronfman's data, as well as clouds within the smallest brightness 
bin from the Georgelin & Georgelin sample. For illustration, the 
thin lines show the locations of the gas spiral arms in our stan- 
dard i/'bar = 20° COBE bulge and disk model without dark halo; 
cf. Section 4. 



and ~ 25° longitude, which is also evident from the Solomon 
et al. (1985) data. 

Table 1 lists a number of tracers that have been used 
to delineate spiral arms. The inferred spiral arm tangents 
coincide with features along the terminal curve in Fig. [jj 
compare also Fig. [l] and Fig. Q. In the inner Galaxy there are 
thus five main arm tangents, of which the Scutum tangent 
is double in a number of tracers. The inner Scutum tangent 
at I ~ 25° is sometimes referred to as the northern 3-kpc 
arm. 

While the main spiral arm tangents on both sides of 
the Galactic center are thus fairly well-determined, it is 
much less certain how to connect the tangents on both sides. 
From the distribution of HII- Regions, Georgelin & Georgelin 
(1976) sketched a four spiral arm pattern ranging from about 
4 kpc galactocentric radius to beyond the solar radius. Al- 
though their original pattern has been slightly modified by 
later work, the principal result of a four-armed spiral struc- 
ture has mostly been endorsed (see the review in Vallee 
1995). 

For illustration, Fig. || also shows the traces in the (I, v) 
diagram of the spiral arms in our standard no-halo model 
(see Section 4) . This model has two pairs of arms emanating 
from the end of the COBE bar, and a four-armed spiral pat- 
tern outside the bar's corotation radius. The model matches 
the observed tangents rather well and illustrates the ways 
in which they can be connected. 



2.5 Dense gas in the Galactic Centre 

Dense regions, where the interstellar medium becomes op- 
tically thick for 21-cm line or 12 CO emission, can still be 
observed using rotational transition lines of 13 CO, CS, and 
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Inner Galaxy spiral arm tangents 
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Table 1. Observed spiral arm tangents compared to model predictions. 



other rare species. In the CO emission line, which probes 
regions with ~ 40 times higher volume densities than the 
12 CO line, an asymmetric parallelogram-like structure be- 
tween ~ +1.5° and ~ —1° in longitude is visible (Bally et al. 
1988). This is nearly coincident with the peak in the termi- 
nal velocity curve and has been associated with the cusped 
xi-orbit by Binney et al. (1991); see also Section 4. In this 
interpretation, part of the asymmetry is accounted for by the 
perspective effects expected for this elongated orbit with a 
bar orientation angle ifha,? = 20°. 

At yet higher densities, the CS line traces massive 
molecular cloud complexes, which are presumably orbit- 
ing on a;2-orbits inside the bar's Inner Lindblad Resonance 
(Stark et al. 1991, Binney et al. 1991). These clouds appear 
to have orbital velocities of < 100 km s~ . 



2.6 Tilt and asymmetry 

We should finally mention two observational facts that are 
not addressed in this paper. Firstly, the gaseous disk be- 
tween the a;2-disk and the 3-kpc-arm is probably tilted out 
of the Galactic plane. Burton & Liszt (1992) give a tilt an- 
gle of about 13° for the HI distribution and show that, by 
combining this with the effects of a varying vertical scale 
height, the observed asymmetry in this region can be ex- 
plained. Heiligman (1987) finds a smaller tilt of ~ 7° for the 
parallelogram. 

Secondly, the molecular gas disk in the Galactic Cen- 
ter is highly asymmetric. Three quarters of the 13 CO and 
CS emission comes from positive longitudes and a differ- 
ent three quarters comes from material at positive velocities 
(Bally et al. 1988). Part of the longitude asymmetry may 
be explained as a perspective effect, and part of both asym- 
metries is caused by the one-sided distribution of the small 
number of giant cloud complexes. Nonetheless it is possible 
that the observed asymmetries signify genuine deviations 
from a triaxially symmetric potential. 



2.7 Solar radius and velocity — comparing real 
and model (l,v) diagrams 

A model calculation results in a velocity field as a func- 
tion of position, with the length-scale set by the distance 
to the Galactic Center assumed in the deprojection of the 
COBE bulge (Binney, Gerhard & Spergel 1997). These au- 
thors took _Ro = 8 kpc, and throughout this paper we will 
use this value in comparing our models to observations. To 
convert model velocity fields into {l,v) diagrams as viewed 
from the LSR, we scale by a constant factor (this gives the 
inferred mass-to-light ratio) and then subtract the line-of- 
sight component of the tangential velocity of the LSR, as- 
suming V = 200kms~\ This value is in the middle of the 
range consistent with various observational data (Sackett 
1997), and is also a reasonable value to use if the Galac- 
tic potential near the Sun is slightly elliptical (Kuijken & 
Tremaine 1994). If the model has a constant circular rota- 
tion curve, i.e., if it includes a dark halo, the LSR velocity 
is part of the model and is scaled together with the gas ve- 
locities. For these models the final scaled LSR tangential 
velocity will be different from Vq = 200 km s -1 and will be 
stated in the text. The radial velocity of the LSR has been 
set to zero throughout this paper. 



3 THE MODELS 

In this Section, we describe in more detail the models that we 
use to study the gas flow in the gravitational potential of the 
Galactic disk and bulge, as inferred from the COBE/DIRBE 
NIR luminosity distribution. In some of these models the 
gravitational field of a dark halo component is added. Self- 
gravit y of the gas and spiral arms are not taken into account 
until § |4.7j . In the following, we first describe our mass model 
as derived from the COBE/DIRBE NIR data (§3.1), then 
the resulting gravitational potential (§3.2) and closed orbit 
structure (§3.3), the assumptions going into the hydrody- 
namical model (§3.4), and finally the main free parameters 
in the model (§3.5). 
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3.1 Mass model from COBE NIR luminosity 
distribution 

The mass distribution in the model is chosen to represent 
the luminous mass distribution as closely as possible. From 
the NIR surface brightness distribution as observed by the 
COBE/DIRBE experiment, Spergel et al. (1996) computed 
dust-corrected NIR maps of the bulge region using a three- 
dimensional dust model. These cleaned maps were depro- 
jected by the non-parametric Lucy-Richardson algorithm of 
Binney & Gerhard (1996) as described in Binney, Gerhard & 
Spergel (1997; BGS). The resulting three-dimensional NIR 
luminosity distributions form the basis of the mass models 
used in this paper. 

The basic assumption that makes the deprojection of 
BGS work is that of eight-fold triaxial symmetry, i.e., the 
luminosity distribution is assumed to be symmetric with re- 
spect to three mutually orthogonal planes. For general ori- 
entation of these planes, a barred bulge will project to a 
surface brightness distribution with a noticeable asymmetry 
signal, due to the perspective effects for an observer at 8 kpc 
distance from the Galactic Center (Blitz & Spergel 1991). 
Vice-versa, if the orientation of the three symmetry planes 
is fixed, the asymmetry signal in the data can be used to 
infer the underlying triaxial density distribution (Binney & 
Gerhard 1996). Because of the assumed symmetry, neither 
spiral structure nor lopsidedness can be recovered by the 
eight-fold algorithm. However, spiral arm features in the 
NIR luminosity may be visible in the residual maps, and 
may appear as symmetrized features in the recovered den- 
sity maps. 

The orientation of the three orthogonal planes is spec- 
ified by two angles. One of these specifies the position of 
the Sun relative to the principal plane of the bulge/bar; this 
angle takes a well-determined (small) value such that the 
Sun is approximately 14 pc above the equatorial plane of 
the inner Galaxy (BGS). The other angle <pbar specifies the 
orientation of the bar major axis in the equatorial plane rel- 
ative to the Sun-Galactic Center line; this angle is not well- 
determined by the projected surface brightness distribution. 
However, for a fixed value of the assumed y>bar, an essentially 
unique model for the recovered 3D luminosity distribution 
results: BGS demonstrated that their deprojection method 
converges to essentially the same solution for different initial 
luminosity distributions used to start the iterations (see also 
Bissantz et al. 1997). As judged from the surface brightness 
residuals, the residual asymmetry map, and the constraint 
that the bar axial ratio should be < 1, admissable values for 
the bar inclination angle i/?bar are in the range of 15° to 35°. 
We will thus investigate gas flow models with <p b ar in this 
range. 

For the favoured (BGS) y>bar = 20°, the deprojected 
luminosity distribution shows an elongated bulge/bar with 
axis ratios 10:6:4 and semi-major axis ~ 2 kpc, surrounded 
by an elliptical disk that extends to ~ 3.5 kpc on the major 
axis and ~ 2 kpc on the minor axis. Outside the bar, the de- 
projected NIR luminosity distribution shows a maximum in 
the emissivity ~ 3 kpc down the minor axis, which appears 
to correspond to the ring-like structure discussed by Kent, 
Dame & Fazio (1991). The nature of this feature is not well 
understood. Possible contributions might come from stars on 
orbits around the Lagrange points (this appears unlikely in 



view of the results of §3.3 and §4.1 below), or from stars on 
xi-orbits outside corotation or on the diamond shaped 1 : 4 
resonant orbits discussed by Athanassoula (1992a) (however, 
the feature is very strong). The most likely interpretation in 
our view, based on §4 below, is that this feature is due to 
incorrectly deprojected (symmetrized) strong spiral arms. If 
this interpretation is correct, then by including these fea- 
tures in our mass model we automatically have a first ap- 
proximation for the contribution of the Galactic spiral arms 
to the gravitational field of the Galaxy. 

In the following, we will model the distribution of lu- 
minous mass in the inner Galaxy by using the deprojected 
DIRBE L-band luminosity distributions for 15° < ipb&r < 
35° and assuming a constant L-band mass-to-light ratio 
M/Ll- The assumption of constant M/Ll may not be en- 
tirely correct if supergiant stars contribute to the NIR lumi- 
nosity in star forming regions in the disk (Rhoads 1998); this 
issue will be investigated and discussed further in §4.1. To 
obtain a mass model for the entire Galaxy, we must extend 
the luminous mass distribution of BGS, by adding a central 
cusp and a model for the outer disk, and (in some cases) 
add a dark halo to the resulting gravitational potential. 

3.1.1 Cusp 

The density distribution of stars near the Galactic Center 
can be modelled as a power law r~ F . From star counts in 
the K-band the exponent p ~ 2.2 ±0.2 for K=6-8 mag stars 
(Catchpole, Whitelock & Glass 1990). The distribution of 
OH/IR stars near the center gives p ~ 2.0 ± 0.2 (Lindqvist, 
Habing & Winnberg 1992). Using radial velocities of the 
OH/IR stars and the assumption of isotropy, Lindqvist et al. 
determined the mass distribution inside ~ 100 pc. The cor- 
responding mass density profile has p ~ 1.5 between ~ 20 pc 
and ~ 100 pc and steepens inside ~ 20 pc. The overall slope 
is approximately that originally found by Becklin & Neuge- 
bauer (1968, p ~ 1.8). 

In the density model obtained from the DIRBE NIR 
data, this central cusp is not recovered because of the lim- 
ited resolution and grid spacing (1.5°) in the dust corrected 
maps of Spergel et al. (1996). The cusp slope of the depro- 
jected model just outside 1.5° moreover depends on that in 
the initial model used to start the Lucy algorithm. To ensure 
that our final density model includes a central cusp similar to 
the observed one, we have therefore adopted the following 
procedure. For the initial model used in the deprojection, 
we have chosen a cusp slope of p = 1.8, in the middle of 
the range found from star counts and mass modelling. This 
gives a power law slope of p = 1.75 in the final deprojected 
density model at around 400 pc. We have then expanded 
the deprojected density in multipoles pi m {r) and have fitted 
power laws to all pio in the radial range 350 — 500 pc. Inside 
350 pc these density multipoles were then replaced by the 
fitted power laws, extrapolating the density inwards. The 
m =fc terms were not changed; they decay to zero at the 
origin. By this modification the mass inside 350 pc is ap- 
proximately doubled. The implied change in mass is small 
compared to the total mass of the bulge and is absorbed in a 
slightly different mass-to-light ratio when scaling the model 
to the observed terminal velocity curve. 
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Figure 3. Contribution of various planar multipoles to the po- 
tential of the standard ^bar = 20° bar model. 
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Figure 4. Rotation curve in the standard </?bar = 20° model 
(solid line), in the model with substracted ring (dashed), and a 
model with constant outer rotation curve (dotted). The rotation 
curves are given for a scaling constant of £ = 1.075. 



3.1.2 Outer disk 

The deprojected luminosity model of BGS gives the den- 
sity in the range < x,y < 5kpc and < z < 1.4 kpc. 
We thus need to use a parametric model for the mass den- 
sity of the Galactic disk outside R — 5 kpc. In this region 
the NIR emission is approximated by the analytic double- 
exponential disk model given by BGS. The least -squares 
fit parameters are Rd = 2.5 kpc for the radial scale-length, 
and zo — 210 pc and z\ = 42 pc for the two vertical scale- 
heights. These parameters are very similar to those obtained 
by Kent, Dame & Fazio (1991) from their SPACELAB data. 
To convert this model for the outer disk luminosity into a 
mass distribution, we have assumed that the disk has the 
same M/Ll as the bulge, because we cannot distinguish be- 
tween the bulge and disk contributions to the NIR emission 
in the deprojected model for the inner Galaxy. 



3.2 Gravitational potential 

From the density model, we can compute the expansion of 
the potential in multipole components $; m (r) and hence the 
decomposition 

$(r, cp) = $o(r) + "M 7 ", <p) cos(2v?) + $ 4 (r, tp) cos(4</?) (2) 

in monopole $rj, quadrupole $2, and octupole $4 terms. 
Higher order terms do not contribute enough to the forces 
to change the gas flow significantly (see Fig. |^), and are 
therefore neglected in the following. The advantage of this 
multipole approximation is that it is economical in terms 
of computer time (no numerical derivatives are needed for 
the force calculations) . The quality of the expansion for the 
forces was tested with an FFT solver. Errors due to the trun- 
cation of the series are typically below 5kms _1 in velocity. 
Because this is less than the sound speed, such errors will 
not significantly affect the gas flow. 

As already mentioned, it is necessary to modify the po- 
tential near the center because of the unresolved central cusp 
in the COBE density distribution. As described above, we 
have replaced the multipole components of the density by 
power law fits inwards of r = 350 pc before computing the 
corresponding multipole expansion of the potential. Since 
the higher order multipoles pio have smaller power law expo- 
nents than the poo term, this implies that the cusp becomes 
gradually spherically symmetric at small radii. We have cho- 
sen this approach because it did not require a specification 
of the shape of the central cusp. Note that without includ- 
ing the modified cusp the gravitational potential would not 
possess x<2 orbits and therefore the resulting gas flow pattern 
would be different. Fig. ^| shows the contribution of the var- 
ious multipoles to the final COBE potential of our standard 
model with tp = 20°. The rotation curve obtained from this 
potential is shown in Fig ^j. 



3.2.1 Dark halo 

If the Galactic disk and bulge are maximal, i.e., if they have 
the maximal mass-to-light ratio compatible with the ter- 
minal velocities measured in the inner Galaxy, then we do 
not require a significant dark halo component in the bar re- 
gion. This may be close to the true situation because even 
with this maximal M/Ll the mass in the disk and bulge fail 
to explain the high optical depth in the bulge microlensing 
data (Udalski et al. 1994, Alcock et al. 1997) by a factor 
> 2 (Bissantz et al. 1997). Thus in our modelling of the bar 
properties we have not included a dark halo component. 

However, for the spiral arms found outside corotation of 
the bar, the dark halo is likely to have some effect. Since we 
only study the gas flow in the galactic plane, the force from 
the dark halo is easy to include without reference to its de- 
tailed density distribution. We simply change the monopole 
moment in the potential directly such that the asymptotic 
rotation curve becomes flat with a specified circular velocity. 
The rotation curve for our flat rotation model is also shown 
in Fig ^. In this model, the halo contribution to the radial 
force at the solar circle is ~ 23%. 
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Figure 5. Left: Effective potential in the standard </?bar = 20° bar model for Qp = 80kms — ^pc - 1 , showing the usual four Lagrangian 
points in the corotation region. Right: For Qp = 55kms kpc . Because the mass peaks in the disk ~ 3kpc down the bar's minor axis 
now contribute significantly to the potential near the increased corotation radius, there are eight Lagrangian points near corotation for 
this pattern speed. 
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Figure 6. Left: Resonance diagram for the standard i/?bar = 20° bar model. Right: Some x\ and X2 orbits in this model for a pattern 
speed of 60kms — ^pc - 1 . The gap between the second and third orbit from outside shows the location of the 1 : 4 resonance. 



3.3 Effective potential, orbits, 
diagram 



and resonance 



The constructed galaxy models have some special properties, 
due to the mass peaks in the disk ~ 3 kpc down the minor 
axis of the bar. In the more common barred galaxy models, 
the effective potential 



$ eff = $ - -n 2 P R 2 



(3) 



in the rotating bar frame contains the usual four Lagrange 
points around corotation and a fifth Lagrange point in the 
centre. In our case, this is true only for larger pattern speeds 
flp, say 80kms _1 kpc _1 , when the corotation region does 
not overlap with the region affected by these (presumably) 
spiral arm features. For lower pattern speeds, in particular 



for Sip ~ 55 — 60kms _1 kpc _1 which we will find below to 
be appropriate for the Milky Way bar, the situation is dif- 
ferent: in this case we obtain four stable and four unstable 
Lagrange points around corotation. The four unstable points 
lie along the principal axes where normally the four usual 
Lagrange points are located, whereas the stable Lagrange 
points lie between these away from the axes (see Fig. ™. For 
yet lower pattern speeds, the number of Lagrange points re- 
duces to four again, but then the two usual saddle points 
have changed into maxima and vice versa. 

We have not studied the orbital structure in this po- 
tential in great detail. However, some of the orbits we have 
found are shown in Fig. [], demonstrating the existence of X\, 
X2, and resonant 1:4 orbits also in this case when there are 
eight Lagrange points near corotation. This is presumably 
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due to the fact that these orbits do not probe the potential 
near corotation. For the orbit nomenclature used here see 
Contopoulos & Papayannopoulos (1980). 

With appropriate scaling, the envelope of the x 1 or- 
bits in our model follows the observed terminal curve in the 
longitude-velocity diagram, and the peak in the curve cor- 
responds approximately to the cusped x\ orbit, as in the 
model of Binney et al. (1991). This will be discussed further 
in §4. 



3.4 Hydrodynamical Models 

For the hydrodynamical models we have used the two- 
dimensional smoothed particles hydrodynamics (SPH) code 
described in Englmaier & Gerhard (1997). The gas flow is 
followed in the gravitational potential of the model galaxy 
as given by the multipole expansion described in §3.2, in a 
frame rotating with a fixed pattern speed Qp . In some later 
simulations we have included the self- grav ity of the spiral 
arms represented by the gas flow (see §4.7 below). 

All models assume point symmetry with respect to the 
centre. This effectively doubles the number of particles and 
leads to a factor of %/2 improvement in linear resolution. We 
have checked that models without this symmetry give the 
same results, as would be expected because the background 
potential dominates the dynamics. 

The hydrodynamic code solves Euler's equation for an 
isothermal gas with an effective sound speed c s : 

^ + (vV)v = - Cs 2 ^-V$. (4) 
at p 



This is based on the results of Cowie (1980) who showed 
that a crude approximation to the ISM dynamics is given 
by an isothermal single fluid description in which, however, 
the isothermal sound speed is not the thermal sound speed, 
but an effective sound speed representing the RMS random 
velocity of the cloud ensemble. 

Using the SPH method to solve Eulers equation has the 
advantage to allow for a spatially adaptive resolution length. 
The smoothing length h, which can be thought of denoting 
the particle size, is adjusted by demanding an approximately 
constant number of particles overlapping a given particle. 
The SPH scheme approximates the fluid quantities by av- 
eraging over neighboring particles and, in order to resolve 
shocks, includes an artificial viscosity. This can be under- 
stood as an additional viscous pressure term which allows 
the pre-shock region to communicate with the post-shock re- 
gion, i.e., to transfer momentum. We have used the standard 
SPH viscosity (Monaghan & Gingold 1983) with standard 
parameters a = 1 and = 2. This SPH method was tested 
for barred galaxy applications by verifying that the prop- 
erties of shocks forming in such models agree with those 
found by Athanassoula (1992b) with a grid-based method. 
See Steinmetz & Miiller (1993) and Englmaier & Gerhard 
(1997) for further details. 

In the low resolution calculations described below we 
have generally used 20000 SPH particles and have taken 
a constant initial surface density inside 7kpc galactocen- 
tric radius. With the assumptions that the gas flow is two- 
dimensional and point -symmetric, these parameters give an 
initial particle separation in the Galactic plane of 62 pc. High 
resolution calculations include up to 100,000 SPH particles 



and may cover a larger range in galactocentric radius to in- 
vestigate the effects of the outer boundary. 

We have experimented with two methods for the initial 
setup of the gas distribution. Method A starts the gas on 
circular orbits in the axisymmetric part $0 of the potential. 
Then the non-axisymmetric part of $ is gradually intro- 
duced within typically one half rotation of the bar. Method 
B places the gas on si-orbits outside and on a;2-orbits in- 
side the cusped a^-orbit. The latter method leads to a more 
quiet start than the former, since the gas configuration is 
already closer to the final equilibrium. Most models shown 
in this paper have been set up with Method A. One model 
was created with a combination of both methods to improve 
resolution around the cusped orbit (see § |l.4| ). 

Different models are usually compared at an evolution- 
ary age of 0.3 Gyr, when the gas flow has become approxi- 
mately quasi-stationary (see §4.2). This corresponds to just 
under three particle rotation periods at a radius of 3kpc. 
The turn-on time of the bar is ~ 0.04 Gyr and is included 
in the quoted evolution age. 

Since the mass-to-light ratio of the model is not known 
a priori, all velocities in the model are known only up to a 
uniform scaling constant f. This implies that also the final 
sound speed is scaled from the value used in the numerical 
calculation: If we know the solution v(r) in one potential $ 
with galaxy mass M, and sound speed c s , then we also know 
the scaled solution £v(r) for the potential £ 2< J>, galaxy mass 
£ 2 M, and sound speed £c s . For the gas we thus effectively 
assume an isothermal equation of state with sound speed 
c s = ^lOkms" 1 . Note that, because c s is a local physical 
quantity, our model can not simply be scaled down to a 
dwarf size galaxy, because then the resulting sound speed 
would be too small and this matters because the gas flow 
pattern depends on this parameter (Englmaier & Gerhard 
1997). 

Below, we fix the scaling constant £ by fitting the ter- 
minal velocity curve of the model to the observed terminal 
velocity curve. To do this, we have to simultaneously assume 
a value for the local standard of rest (LSR) circular motion. 
These two parameters compensate to some extent, but we 
generally find that the fit to the terminal curve is more sen- 
sitive to the assumed LSR motion than to the value of the 
scaling constant. In most of the rest of the paper we have 
therefore fixed the LSR velocity to 200 km s -1 . We work in 
units of kpc, Gyr, and Mq. If the deprojected COBE den- 
sity distribution is assumed to be in units of 3 1O 8 M0 kpc -3 , 
£ is found to have typical values of 1.075 to 1.12. 

3.5 Summary of model parameters and discussion 
of assumptions 

Our models have a small number of free parameters; these 
and the subset which are varied in this paper are listed here. 
This Section also contains a brief summary and discussion 
of the main assumptions. 



3. 5. 1 Bar parameters 

• Probably the most important parameter is the bar's coro- 
tation radius R c or pattern speed f2p, which will set the 
location of the resonance radii and spiral arms and shocks 
in the gas flow. 
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• The second important bar parameter is its orientation an- 
gle ip with respect to the Sun-Galactic Center line, which 
affects the appearance of the gas flow as viewed from the 
Sun. 

The shape and radial density distribution in the bar re- 
gion are constrained by the observed NIR light distribution. 
However, their detailed form is dependent on the assumption 
of 8-fold symmetry, which is likely to be a good assumption 
in the central bulge region, but might be too strong in the 
outer bar regions, where a possible spiral density wave might 
affect the dynamics. It is also possible that an overall m = 1 
perturbation is needed to explain the observed asymmetries, 
such as in the distribution of giant cloud complexes in the 
Galactic center, or the fact that the 3-kpc-arm appears to be 
much stronger than its counterarm. Nevertheless, it is im- 
portant to find out how far we can go without these asymme- 
tries. In any case, the bar should have the strongest impact 
on the dynamics. 



3.5.2 Mass model 

• The only additional parameter in the luminous mass model 
is the scaling constant £ which relates NIR luminosity and 
mass. For each pair of values of the previous two parameters 
and at fixed LSR rotation velocity this is determined from 
the Galactic terminal velocity curve, assuming that this is 
dominated by the luminous mass in the central few kpc. 

This contains the additional assumption that all compo- 
nents have the same constant NIR mass-to-light ratio. This 
appears to be a reasonably good assumption on the basis 
of the fact that optical-NIR colours of bulges and disks in 
external galaxies are very similar (Peletier & Balcells 1996). 
It is unlikely to be strictly correct, however, because the 
bulge and disk stars will not all have formed at the same 
time. To relax this assumption requires additional assump- 
tions about the distinction between disk and bulge stellar 
luminosity. This is presently impractical. 

• Depending on the LSR rotational velocity, a dark halo is 
required beyond R ~ 5 kpc. Thus we need to specify the 
asymptotic circular velocity of the halo. Here we consider 
only two cases, one without halo, the other with asymptotic 
halo velocity of vo = 200 kms -1 . This is in the middle of the 
observed range of 180 - 220 kms" 1 (Sackett 1997). 



3.5.3 LSR motion and position 

• We assume throughout this paper that the distance of the 
LSR to the galactic centre is Ro — 8 kpc (see the review by 
Sackett 1997, but also the recent study by Oiling & Merri- 
field 1998, who argue for a somewhat smaller Ro). 

To compare model velocity fields with observations, we 
have to know not only the position of the Sun but also its 
motion. The peculiar motion of the Sun relative to the LSR 
is often already corrected for in the published data. 

• The remaining free parameter is the LSR rotational ve- 
locity around the galactic center, which lies in the range 
between 180 - 220 kms -1 (Sackett 1997). We will again use 
Vo = 200 kms -1 , consistent with the above. 



3.5.4 Gas model 

We use a crude approximation to the ISM dynamics, that of 
an isothermal single fluid (Cowie 1980). The effective sound 
speed c s is the cloud-cloud velocity dispersion; this varies 
from ~ 6 kms -1 in the solar neighbourhood to ~ 25 kms - 
in the Galactic Center gas disk. We have considered models 
with globally constant value of c s between 5 — 30 kms -1 and 
have not found any interesting effects. Only at the largest 
values do the spiral arm shocks become very weak. 

Consistent with the assumption of eight-fold symmetry 
for the mass distribution, we have assumed the gas flow to be 
point symmetric with respect to the origin. For gas flows in 
the eight-fold symmetric potential and without self-gravity 
there are no significant differences to the case when the gas 
model is run without symmetry constraint. 



4 RESULTS 

4.1 Gas flow morphology implied by the COBE 
luminosity /mass distribution 

We begin by describing the morphology of gas flows in the 
COBE-constrained potentials. For our starting model we 
take the deprojected eight-fold symmetric luminosity distri- 
bution obtained from the cleaned COBE L-band data, for a 
bar angle iy9bar = 20° as favoured by BGS. This model, with 
constant mass-to-light ratio and no additional dark halo, 
will be referred to as the (standard) ipb aT = 20° COBE bar. 

In our first simulation this bar model is assumed to ro- 
tate at a constant pattern speed il p such that corotation 
is at a galactocentric radius of approximately 3.1 kpc. With 
the value for the mass-to-light ratio Tl as determined by 
Bissantz et al. (1997) from fitting the observed terminal ve- 
locities, this gives Q p = 60 kms -1 /kpc. 

For the gas model we take a constant initial surface den- 
sity on circular orbits, represented by 20000 SPH particles, 
and an effective isothermal sound speed of c s = 10 kms -1 . 
The gas is relaxed in the bar potential as described in §3.4 
(Method A) and the initial particle separation in the Galac- 
tic plane is 62 pc. 

Fig. shows the morphology of the gas flow in this 
model. Inside corotation (R c = 3.1 kpc) two arms arise near 
each end of the bar. The dust lane shocks further in are 
barely resolved in this model. The structure into the coro- 
tation region is complicated. Outside corotation four spiral 
arms are seen. This four-armed structure is characteristic 
of all gas flow models that we have computed in the un- 
modified COBE potentials. A four-armed spiral pattern is 
considered by many papers the most likely interpretation of 
the observational material regarding the various spiral arm 
tracers and the five main spiral arm tangent points inside 
the solar circle (cf. Figs, [l], ^| Table [lj), see the review and 
references in Vallee (1995). 

Two of the four arms emanate approximately near the 
major axis of the bar potential, two originate from near the 
minor axis. The additional pair of arms compared to more 
standard configurations is caused by the octupole term in 
the potential; in a model where this term is removed, the 
resulting gas flow has only two arms outside corotation. 

Fig. M shows the gas flow in a model in which the density 
multipoles with m/0 were set to zero outside 3 kpc before 
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Figure 7. Gas flow in the i^bar = 20° COBE bar with corotation 
at R c ~ 3.1 kpc. In this and subsequent figures, the long axis 
of the bar lies along the x-axis, and the location of the Sun is 
at x = — 7.5kpc, y = — 2.7kpc, 20° away from the bar's major 
axis. The simulation has TV = 20000 SPH particles with point 
symmetry built in. The initial gas disk extends to i? max = 7 kpc. 
The multipolc expansion of the potential includes all significant 
terms (up to I = 6, m = 4). 
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Figure 8. Gas flow in the same model as in Fig. |7J, but with the 
density multipoles with m ^ set to zero outside 3 kpc. Since the 
I = terms are unchanged, the circular rotation curve remains the 
same. The gas model extends to 10 kpc and has 20'000 particles. 



computing the potential. This modification leaves the circu- 
lar rotation curve of the model unchanged. All structure in 
the resulting gas flow is now driven by the rotating bar inside 
corotation, whose quadrupole moment outside R c is weak. 
The figure shows that, correspondingly, only two weak spiral 
arms now form in the disk outside corotation. In the (l,v) 
diagram, these appear as tangents at longitudes I ~ —50° 
and I ~ 50° . However, there are no arms in this model which 
would show along the tangent directions I — ±30°; at best 
there are slight density enhancements in these parts of the 
disk. However, the abundance of warm CO clouds found 
near I — 25 — 30° by Solomon, Sanders & Rivolo (1985; cf. 
Fig. ^) indicates that a spiral arm shock must be present in 
this region. Thus the model underlying Fig. |^, in which all 
structure in the gas disk outside 3 kpc is driven by only the 
rotating bar in the inner Galaxy, cannot be correct. 

Both the quadrupole and octupole terms of the poten- 
tial outside ~ 3 kpc are dominated by the strong luminosity- 
mass peaks about 3 kpc down the minor axis of the COBE 
bar. From comparing Figs. |^ and ^ we thus conclude that, 
in order to generate a spiral arm pattern in the range 
R — 3 — 8 kpc in agreement with observations these peaks 
in the NIR luminosity must have significant mass. In other 
words, the NIR mass-to-light ratio in this region cannot be 
much smaller than the overall value in the bulge and disk. 
This result is in agreement with a recent study by Rhoads 
(1998) who finds that in external galaxies the local contri- 
bution of young supergiant stars to the NIR flux can be of 
order ~ 33% but does not dominate the old stellar popula- 
tion. 

The observed luminosity peaks on the minor axis of the 
deprojected COBE bar coincide with dense concentrations of 
gas particles near the heads of the two strongest spiral arms 
in our gas model (Fig. Q). Since the gas arms will generally 
be accompanied by stellar spiral arms, this suggests that the 
most likely interpretation of the minor axis peaks in the BGS 
model is in terms of incorrectly deprojected spiral arms. 

Spiral arms are generally the sites of the most vigorous 
star formation in disk galaxies. Also in the Milky Way, ob- 
servations of far infrared emission show that most of the star 
formation presently occurs in the molecular ring (Bronfman 
1992). Since we have found from dynamics that even in this 
region the associated young supergiants do not dominate the 
NIR light, this implies that over most of the Galactic disk 
the assumption of constant NIR mass-to-light ratio for the 
old stars is justified. 

4.2 Time evolution 

How stationary is the morphological structure in these gas 
flows? To address this question, we show in Fig. ^ the 
time evolution of a typical model ((£>bar = 20°, fip = 
55 kms _1 kpc _1 ). In this and other simulations the non- 
axisymmetric part of the gravitational potential was gradu- 
ally turned on within about one half of a bar rotation period 
(~ 0.04 Gyr). The gas flow, which is initially on circular 
orbits, then takes some time to adjust to the new poten- 
tial. It reaches a quasi-stationary pattern by about time 
t = 0.3 Gyr. This flow is shown in the top left panel of 
Fig. After t = 0.3 Gyr the variations in the gas flow are 
small: about 5kms _1 in the velocities. Also the sharpness 
of the arms inside corotation varies slightly. In this quasi- 
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Figure 9. Gas particle distribution in the standard <£bar = 20° COBE bar potential, for a corotation radius R c ~ 3.4 kpc. The frames 
show snapshots at t = 0.3 Gyr, 1.0 Gyr, 2.0 Gyr, 3.0 Gyr (top left, top right, bottom left, bottom right). The bar was gradually turned on 
between t = and t = 0.04 Gyr. The particle rotation period at 3 kpc galactocentric radius is ~ 0. 1 Gyr. The most significant evolutionary 
effect is a loss of resolution by about a factor of two in linear distance between the first and the last frame, due to substantial mass 
inflow. This results in fuzzier spiral arms at the end of the simulation, which appear to terminate earlier. Also compare Fig. h 01 



stationary flow material continuously streams inwards: Gas 
parcels that reach a shock dissipate their kinetic energy per- 
pendicular to the shock. Subsequently they move inwards 
along the shock. 

As Fig. |^ shows, the inward gas inflow causes a slow 
evolution without much changing the morpophology of the 
gas flow. However, the mass accumulating on the central 
disk of a^-orbits in the course of this process is consider- 
able. In fact, to continue the simulation we have found it 
necessary to constantly remove particles from the a;2~disk. 
In doing this we have simultaneously increased the particle 
mass in this region in such a way as to keep the surface den- 
sity unchanged. Therefore effectively we have only limited 



the resolution in this region from increasing ever further, 
without rearranging or changing mass. The gas inflow leads 
to a loss of resolution in the outer disk. From t — 0.3 Gyr 
to t = 3 Gyr, the surface density of particles in outer disk 
of the model shown in Fig. ^ decreases by about a factor 
of four, i.e., the linear resolution by about a factor of two. 
The spiral arms therefore become more difficult to see; in 
particular, the starting points and end points of some of the 
arms appear to shift slightly. 

The most rapid evolution occurs in the vicinity of the 
cusped orbit. Already by time t = 0.3 Gyr, the gas disk near 
this orbit has been strongly depleted. Because the shear in 
the velocity field in the vicinity of the cusped orbit is very 
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Figure 10. Longitude-velocity (l,v) diagrams corresponding to 
the gas particle distributions in Fig. H at t = 0.3 Gyr (top) 
and 3 Gyr (bottom). In constructing these we have assumed 
Ro = 8kpc and vo = 200kms — 1 . The inner disk on X2-orbits, 
the terminal velocity curve, and the spiral arm traces are appar- 
ent. 



strong, particles that reach the cusped orbit shock move 
to the center quickly along the shock ridges and then fall 
onto the £2-disk (Englmaier & Gerhard 1997). In the low- 
density region between the cusped orbit and the a;2-disk, the 
smoothing radius of the SPH particles is large compared to 
the velocity gradient scale. It is possible that the result- 
ing large effective viscosity accelerates the depletion of gas 
on the cusped orbit and in its vicinity. A similar effect has 
been seen by Jenkins & Binney (1994) in their sticky parti- 
cle simulations. See § 4.4 for an improved model with more 



resolution in the cusped orbit region. 



4.3 Model (l,v)-diagrams 

The non-axisymmetric structures seen in Figs. ^ etc. lead 
to perturbations of the gas flow velocities away from cir- 
cular orbit velocities. These can be conveniently displayed 
in an (l,v) diagram like those often used for representing 
Galactic radio observations. In fact, to constrain the Galac- 
tic spiral arm morphology from comparisons of our models 
with Galactic radio observations, we really only have (l,v) 
diagrams! Fig. |lo| shows (I, v) diagrams obtained from the 
gas distributions in the first and last panel of Fig. ^, at 
t = 0.3 Gyr and t = 3 Gyr, respectively, for an assumed dis- 
tance of the Sun to the galactic center of 8.0 kpc and LSR 
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Figure 11. Top: Schematic representation of the spiral arms in 
the gas flow model depicted in the top left frame of Fig. ^. 
Bottom: The same spiral arms in the (I, v) diagram. 



rotation velocity vq — 200 km s -1 . The bright ridge rising 
steeply from the center in these diagrams is caused by the 
dense disk of gas on a^-orbits, visible in the very center of 
the flow in Fig. |^. The more irregularly-shaped ridges are 
the traces of spiral arms in the (I, v) diagram. Also well vis- 
ible in Fig. [id are the terminal velocity curves. 

Looking at Figs. ^ and [lo] shows that the relation be- 
tween morphological structures in the gas disk and corre- 
sponding structures in the (l,v) diagram is somewhat non- 
intuitive. In order to gain a better understanding of this 
relation, we have constructed a schematic representation 
of the arm structures of the model in the top left panel 
of Fig. ^ in both coordinate planes. The upper panel in 
Fig. |ll] shows schematically the location of the gaseous spi- 
ral arms, the cusped orbit shocks (dust lanes), and the xi- 
disk in the (a;, j/)-plane. In this diagram, the Sun is located 
at x — —7.5 kpc, y — — 2.7 kpc, i.e. at Ro — 8 kpc and 
<^bar = 20°. The lower panel shows the corresponding fea- 
tures in the (l,v) diagram as observed from this LSR posi- 
tion, with the same line styles to facilitate cross identifica- 
tion. 
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We see that whenever a spiral arm crosses a line-of- 
sight from the Sun twice, it appears as a part of a loop in 
the (l,v) diagram. This is the case, e.g., for the two outer 
spiral arms seen nearly end-on (thin full lines in Fig. |l"lj ), 
and for the innermost pair of arms driven by the bar (thick 
and thin dashed lines). The equivalent to the 3kpc arm (see 
below) and the corresponding counterarm on the far side 
of the galaxy are parts of a second pair of arms driven by 
the bar; these cross the relevant lines-of-sight to the Sun 
only once (thick full lines and thick dotted lines in Fig. [y]) , 
respectively). The same is true for the outer pair of spiral 
arms seen nearly broad-on as viewed from the Sun (dash- 
dotted and small dotted lines). 

It is clear from inspection of Figs. 0and|ll| and a com- 
parison with the corresponding observational data (see the 
figures reproduced in Section 2, and the diagrams in the pa- 
pers cited there) that already the initial COBE-constrained 
model gas flow of Fig. ^| resembles the Milky Way gas dis- 
tribution in several respects: 

(i) The number of arm features in the longitude range 
[—60°, 60°] and their spacing in longitude is approximately 
correct (compare Table |l]). 

(ii) The model contains an arm which passes through 
the I = 0-axis at negative velocity (~ — 30kms _1 ) and 
merges into the southern terminal velocity curve at negative 
I. Qualitatively, this is similar to the well-known 3kpc-arm, 
although this crosses the I = 0-axis at ~ — 50kms _1 and 
extends to larger longitudes. 

(iii) The positions and velocities of gas particles in the 
22-disk are similar to those observed for the giant Galactic 
Center molecular clouds on the I > side in the CS line 
(Bally et al. 1987, 1988; Binney et al. 1991). 

(iv) The terminal velocity curve slopes upwards towards 
large velocities near I — 0, although not as much as would 
be expected for gas on and just outside the cusped orbit in 
the COBE potential (see § 4.4), and not as much as seen 
in the HI and CO data. In the following subsections we 
compare the COBE models more quantitatively with the 
observational data, and attempt to constrain their main pa- 
rameters. 



4.4 The terminal velocity curve 

We will now determine the mass normalisation of the models 
by fitting their predicted terminal velocity curves to obser- 
vations. Model terminal curves are computed by searching 
for the maximal radial velocity along each line of sight, as 
seen from the position of the Sun 8 kpc from the center, and 
then substracting the projected component of the LSR mo- 
tion. For the latter we have assumed that the LSR motion is 
along a circular orbit with wlsr = 200kms _1 and no radial 
velocity component. By an eye-ball fit of these model curves 
to the observed terminal velocit y cu rve we then obtain the 
proper scaling constant £ (see 



3.4) by which all velocities 



in the model have to be multiplied to obtain their Galactic 
values. Both the mass and the potential then scale with £ 2 . 

Figure ^ shows several model terminal velocity curves 
obtained in this way, and compares them with the north- 
ern and southern Galactic terminal velocities as determined 
from HI and 12 CO observations from a number of sources. 
Although the error bars for all measured terminal velocities 
are small, the data show some scatter due to differences in 



i — i — i — I — i — i — i — 3c — I — i — i — i — i — I — i — i — r 




Figure 13. Top: The upper envelope of the a'l-orbits in the (I, v)- 
plane traces out the observed inner terminal velocity curve. The 
figure shows xi-orbits computed in the same gravitational poten- 
tial as that used for the models shown in Fig. |l^; each xi -orbit 
traces a parallelogram-like curve in this plot. The innermost x\- 
orbit shown is the cusped orbit which peaks at I = 2°. Inside 
this orbit, the gas switches to 22-orbits which reach velocities 
of 120 kms~ 1 in this potential. The observed northern terminal 
velocities are shown for comparison (same sources as in Fig. |l2"| ). 

Bottom: Particle distribution in the same part of the (I, v)- 
plane, for the high— resolution model shown by the thick solid line 
in Fig. [l2| The cusped orbit and the innermost few x\ -orbits are 
not occupied even at this resolution. The parallelogram structure 
in this plot therefore reaches inwards and upwards to only I ~ 3° 
and v ~ 220 kms -1 , compared to I ~ 2° and v ~ 270kms — 1 
for the orbits, and / ~ 2° and v ~ 260 kms -1 for the observed 
terminal velocities. 



angular resolution and sensitivity. The model terminal ve- 
locity curves Fig. [l^ are for models of different spatial res- 
olution and simulated radial extent, to illustrate the effects 
of these parameters. The dash-dotted line is for a model 
extending to galactocentric radius 9 kpc; this demonstrates 
that for wlsr > 200 km s^ 1 the observed terminal curve be- 
yond about ±(40 — 50)° requires a dark halo component in 
the Milky Way. 

After including a halo component in the model (by mod- 
ifying the monopole component of the potential such that 
the rotation velocity becomes constant outside 4.5 kpc; see 
Fig. ^) , both the northern and southern terminal curves are 
much better reproduced (thick solid line in Fig. 0). The 
rotation curve of this halo model is shown in Fig. Q the 
simulated radial range is 12 kpc. The contribution of the 
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Figure 12. Northern and southern Galactic terminal velocity curves compared with model predictions. Observational data from sources 
as follows. Filled squares (1 = — 10°): HI data from Fig. 1 of Burton & Liszt (1993). Empty squares (1 = — 20°): unpublished 
140ft single dish HI data, kindly provided by Dr. B. Burton. Open circles: HI terminal velocities from Fich et al. (1989), based on data 
from Westcrhout (1957). Diamonds, with tiny error bars: northern 12 CO terminal velocities from Clemens (1985). Without error bars: 
southern 12 CO data from Alvarez et al. (1990). The data are relative to the LSR, mostly as corrected by the respective authors. The 
Clemens (1985) data have been corrected for internal dispersion, and the velocities in Burton & Liszt (1993) are relative to the Sun and 
have been corrected for LSR motion. 

The model terminal velocity curves are from: a typical low resolution model with no halo and declining rotation curve (dash-dotted); 
the same model at higher resolution but with less radial extent (dotted); a low resolution model with fiat rotation curve (thin solid); and 
our best model with high resolution (see text) and flat rotation curve (thick solid line). All models assume a bar angle i^bar = 20° and 
have corotation at R c ~ 3.4 kpc. 
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dark halo inside the solar radius is fairly small (~ 23% in 
the radial force at 8kpc), somewhat less even than in Kent's 
(1992) maximum disk model. 

In this model, there remain two main regions of dis- 
crepancy with the observed terminal velocity curve. First, 
the model terminal velocities are too low at and just out- 
side the peak at Z ~ 2°. This is strongly influenced by and 
probably due to resolution effects, as discussed below. Sec- 
ond, there is a larger mismatch around —20°. This is likely 
caused by our mass model not being correct in the vicinity 
of the NIR lumps ~ 3kpc down the minor axis of the bar. 
Lines-of-sight at around —20° cross one of these lumps as 
well as the end of the 3 kpc arm and the head of one of the 
spiral arms outside corotation (see Fig. pd| ). The eight-fold 
symmetric deprojection of BGS is therefore likely to give in- 
correct results in this region. Smaller systematic deviations 
in the terminal velocities are visible around Z = 30 — 50° and 
I — — (50 — 70)°, although there, and everywhere else, the 
differences between model and observations are now of the 
order of the scatter between the various observational data 
and of the order expected from perturbations in the disk. 
Given the uncertainties the overall agreement is surprisingly 
good. This suggests that the basic underlying assumption, 
that in the inner Galaxy the NIR light traces the mass, is 
mostly correct. 

Fitting the terminal velocity curve to both sides we ob- 
tain a scaling constant of £ = 1.12. This is slightly larger 
than the value obtained in our first attempts to fit the mod- 
els to the observations (£ = 1.075, Bissantz et al. 1997), in 
which we only considered the northern rotation curve and 
ignored the data beyond I — 48°. It is worthwhile pointing 
out that the derived value of £ is only weakly dependent on 
the assumed LSR tangential velocity: In order of magnitude, 
a 10% change in the LSR velocity leads to a 1% difference 
in f. With an improved dust model and deprojection of the 
outer disk we could therefore attempt to determine Vo from 
these models. 

As is clearly visible in Fig. |l2| the peak in the observed 
terminal curve at 2° is not well reproduced by our lower- 
resolution models. However, Fig. [llfl shows that it is nicely 
approximated by the envelope of the a;i-orbit family when 
all orbital velocities are scaled by the same value of f. At 
early times in the model evolution, when the gas flow is 
not yet stationary, the peak is also reproduced in the hy- 
drodynamic gas model, but thereafter the region around 
the cusped orbit is depopulated (see also Jenkins & Bin- 
ney 1994). We attribute this to the artifical viscosity in the 
SPH method, which smears out the velocity gradient over 
two smoothing lengths, and to the method used for setting 
up the gas simulation. 

We can estimate the magnitude of the effect as follows. 
Near the cusped orbit, which sets the maximum velocity 
along the terminal velocity curve, the particle smoothing 
length h in the low resolution model is large, about ~ 100 pc, 
because the gas density in this region is small. The full x\— 
orbit velocity on the terminal curve can only be reached 
about two smoothing lengths away from the cusped orbit, 
where it is no longer affected by the more slowly moving 
gas on X2-orbits further in. The longitudinal angle corre- 
sponding to about 2h at the distance of the galactic centre 
is about 1.4°. Thus, the peak in the gas dynamical terminal 



curve should be found at / > 3.4° when the orbits peak at 
2°, showing how sensitive the peak location is to resolution. 

To test this explanation, we have run a bi-symmetric 
model with 100, 000 particles, resulting in about 2.2 times as 
much spatial resolution as in the low-resolution models. Just 
on the basis of this higher resolution, the terminal velocity 
peak then moves from about 5° to 4° (at t — 0.3 Gyr). We 
then further increased the resolution by the following proce- 
dure. The gas inside the outermost x± orbit shown in Fig. ^ 
was removed, and set up again on nested closed x\ and X2 
orbits, while keeping particles outside this region unchanged. 
Evolving this modified gas distribution for a further 0.3 Gyr, 
we obtained our final high resolution model. This is shown 
by the thick solid line in Fig. which peaks at about 3° 
and vt — 235 km s" 1 . Compared to the original 20000 par- 
ticle model (thin solid line in Fig. [li]) the mismatch at the 
peak has been reduced by about a factor of two in scale and 
by two thirds in the peak velocity. 

Although this analysis was inspired by a technical prob- 
lem, there is an observable implication of it as well. Since 
we may interpret the hydrodynamical model in terms of gas 
clouds having a mean free path length of order the smooth- 
ing length, we may restate the result in the following way: 
a loss of resolution occurs when the cloud mean free path 
is significant compared to the gradient in the true velocity 
field. Applied to the inner Galaxy, our result then indicates 
that the clouds near the cusped orbit peak in the termi- 
nal velocity curve must have short mean free paths, i.e., be 
described well in a fluid approximation. 

Apart from resolution effects, the precise position of the 
peak in the terminal velocity curve also depends critically 
on the location of the ILR and hence on the mass model in 
the central few 100 pc. In this region, the deprojected COBE 
model suffers from a lack of resolution and our added nuclear 
component has uncertainties as well. We therefore believe 
that with improved data and further work the remaining 
discrepancies in this region will be fixed. 

4.5 Pattern Speed and Orientation of the 
Galactic Bar 

There are two observations which constrain the value of the 
pattern speed rather tightly. First, there is the 3-kpc-arm, 
a feature which exhibits non-circular motions of at least 
50kms _1 . In our models, we find that only the arms in- 
side the bar's corotation radius are associated with strong 
non-circular motions, so such an arm has to be driven by 
the bar. From observations and models of barred galaxies 
we also know that strong spiral arms associated with both 
ends of a bar are common. Therefore we conclude that the 
3-kpc-arm must lie inside the bar's corotation radius. 

Second, a lower limit to the corotation radius is given 
by the inner edge of the molecular ring. If the molecular 
ring were indeed a ring such as induced by a resonance, it 
would be located near the outer Lindblad resonance (e.g., 
Schwarz 1981). On the other hand, if it is actually made of 
several spiral arms (Dame 1993, Vallee 1995, this paper), 
then the small observed non-circular velocities along these 
spiral arms also show that these arms must be outside the 
bar's corotation radius. Solomon et al. (1985) find from the 
distribution of hot, presumely shocked cloud cores, that the 
inner edge of the molecular ring is at R — 4 kpc. The total 
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Figure 14. Longitude-velocity (I, v) diagrams for the gas flows in COBE bar potentials with different pattern speeds and bar orientation 
angles. The column of (l,v) diagrams on the left shows the influence of the pattern speed on the gas flow in the standard ^bar = 20° 
COBE bar potential. These frames show models with corotation radii R c = 3.1 kpc (top) R c = 3.4kpc (middle) to R c = 4.0kpc (bottom). 
The right column shows gas flows in COBE bars deprojected for different bar orientation angles, i^bar = 15° (top), i/?bar = 25° (middle), 
Vbar = 30° (bottom), all for corotation at R c = 3.4 kpc. 



surface density of neutral gas also drops dramatically inside 
of 4 kpc (Dame 1993). From the IRT 2.4 jj,m photometry of 
the galactic disk Kent et al. (1991) concluded that there 
is a ring, or spiral arm, at about R = 3.7 kpc. Therefore 
we conclude that the bar's corotation radius is inside R = 
4 kpc. 

An independent argument for corotation falling some- 
where between 3 kpc and 4 kpc comes from the fact that the 
deprojected COBE bar appears to end somewhere between 
3 kpc and 3.5 kpc (BGS). From both N-body simulations 
and direct and indirect observational evidence the corota- 



tion radius is usually found at between 1.0 and 1.2 times 
the bar length (Sellwood & Wilkinson 1993; Merrifield & 
Kuijken 1995; Athanassoula 1992b). In our models, a cora- 
tion radius between 3 kpc and 4 kpc corresponds to a pattern 
speed of ~ 50 — 60kms _1 kpc _1 . 

We have run gas dynamical simulations with corota- 
tion at 4.0, 3.4, and 3.1 kpc, to determine from observations 
which of these values is most nearly appropriate. For the 
comparison with observations it is important to notice that 
several other parameters enter here, most importantly, the 
orientation angle of the bar, the uncertain contribution of 
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the dark halo to the outer rotation curve and hence terminal 
velocities, and the LSR velocity. 

We first fix the bar orientation angle at y>bar = 20°, 
but will vary this parameter later. The chosen value of vJbar 
is in the range allowed by the NIR photometry (BGS), it is 
favoured by the clump giant star distribution as analyzed by 
Stanek et al. (f997) and by the gas kinematical analysis of 
the molecular parallelogram by Binney et al. (1991), and it 
meets the preference for an end-on bar in the interpretation 
of the microlensing experiments. 

In the last section we found that the Galactic terminal 
velocity curve for \l\ < 45° is well-reproduced by the gas 
flow in the maximum NIR disk model with constant mass- 
to-light ratio. Moreover, even this maximum disk model fails 
by a factor of > 2 in explaining the high microlensing opti- 
cal depth towards the bulge (Bissantz et al. 1997), making 
it very difficult to further reduce the mass in the interven- 
ing disk and bulge. We can therefore confidently assume a 
maximum disk model in the following and, to separate the 
determination of the bar and halo parameters, we restrict 
the comparison with observations to longitudes |Z| < 45°. 

Finally, we set the LSR rotation velocity to Vo = 
200 km s -1 , in the middle of the observed range (§2.7). A 
10% difference in this parameter is not very important for 
the comparison with the inner Galaxy gas velocities. 

Thus we begin by considering a sequence of models 
with varying corotation radius R c and the other parameters 
fixed as just described. For three models with R c — 4.0 kpc, 
3.4 kpc, and 3.1 kpc we have plotted (l,v) diagrams and have 
determined the scaling constant £ for each simulation by fit- 
ting to both terminal curves. The final scaled (I, v) diagrams 
are shown in the left column of Fig. |l4]. For the scaling con- 
stant we obtain £ = 1.13, 1.12, and 1.09 for the 4.0, 3.4, 
and 3.1 kpc models. The correctly scaled pattern speeds are 
then 59, 57, and 61 kms _1 kpc _1 . This means that by chang- 
ing the corotation radius, we effectively change the mass of 
the model galaxy, while keeping the pattern speed almost 
constant at about 60kms~ 1 kpc _1 . 

At small absolute longitudes, \l\ < 10°, these low- 
resolution models do not have enough particles to resolve 
the true gas flow, and furthermore there are no published 
terminal velocities in this region on the southern side. Thus 
for now we ignore data near the peak of the terminal velocity 
curve. This leaves a range of ±Z = 10 — 45° within which we 
compare these no-halo model gas flows with the northern 
and southern terminal velocity curves and with the various 
spiral arm features shown in Figs, [l] ^| 

On the northern side (I > 0), we try to match the model 
to the pronounced spiral arm at about +30°, which is best 
visible in the warm CO clouds (Solomon et al. 1985) and in 
the distribution of HH-regions. Moreover, the (l,v)-diagrams 
of CO and HH-regions shows that the +30° arm is double 
(Fig. |^). In the model shown schematically in Fig. |ll], there 
are actually three arms near I ~ 30°, two of which overlap, 
while the third, the northern 3-kpc-arm (thick dotted line), 
runs almost parallel to the first two. There is also a wiggle in 
the terminal curve at about +10°, which is probably caused 
by a spiral arm similar to the northern, secondary inner 
spiral arm in the model (thin dashed line in Fig. |llj). The 
southern terminal curve is more distorted by spiral arms 
than the northern curve. A pronounced feature is the —30° 
arm in the molecular ring, as well as the well-known 3-kpc- 



arm which continues on from a non-circular velocity ridge 
beginning at I ~ 10° and v = 0. 

Fig. u4 (left column) shows that the gas flows in all three 
cases are similar; nonetheless small differences in the spiral 
arm locations help to show that the 3.4 kpc case is closest 
to the real Galaxy. The 30° spiral arm tangent is best re- 
produced in the 3.4kpc-model (middle panel in left column 
of Fig. ^). In the top panel, it is not double as observed, 
and in the bottom panel the tangent moves out to ~ 40°. 
The 50° spiral arm tangent is reasonably well reproduced in 
the top two panels, but is absent or very weak in the bot- 
tom panel, but this arm may not be a reliable indicator in 
the absence of a halo. In the south, we observe that large 
corotation radii move the arm at —30° outwards. The —30° 
spiral arm tangent is in about the correct location in the top 
two panels, but too far out in the bottom panel. The mis- 
match of the terminal curve at about —20° becomes larger 
with decreasing corotation radius as well. On the other hand, 
the —50° spiral arm tangent is not present for the largest 
pattern speed (top panel), while it is adequately present in 
the lower two panels. The 3-kpc-arm equivalent is present 
in all three models, but its detailed locus in the (l,v) dia- 
gram is not correct in any of the models: it has either too 
low non-circular velocities at / = (particularly in the bot- 
tom panel), or it does not extend to large enough negative 
longitudes (particularly in the top panel). For corotation at 
3.4 kpc, the 3-kpc-arm ends in the model exactly at 3 kpc, as 
seen from the assumed solar position. There are particles in 
the low intensity forbidden velocity region bounded by the 
line from (l,v) = (10,0) — (0,-50) in all three models, but 
the resolution does not suffice to prefer one model over the 
other. Finally, an equivalent to the 155kms _1 arm at small 
negative I is present only for the lowest pattern speed. 

From an unweighted average of these comparisons to 
various observational landmarks we conclude that the coro- 
tation radius of the Galactic bar is most likely at about 
3.4 kpc. However, since none of the above models is exactly 
right yet, this value could well change by 10% when other 
effects are taken into acc oun t, such as bar orientation (dis- 
cussed next), dark halo (§4.6), or self-gravitating spiral arms 



(§4.7) 



With the corotation radius fixed at 3.4 kpc, we can at- 
tempt to find an optimal orientation angle for the bar. For 
this purpose, a series of deprojected bar models have been 
made as in BGS, for bar orientation angles i^bar = 10, 15, 20, 
25, and 30°. Smaller or larger angles are not consistent with 
the asymmetry pattern of the observed NIR distribution or 
result in unphysical bar shapes, see BGS and Bissantz et al. 
(1997). With the bar orientation the shape and radial extent 
of the deprojected bar change. 

Gas flow simulations with R c — 3.4 kpc were made for 
each of these cases, and (I, v) diagrams for ip^av = 15, 25, and 
30° are shown in the right column of Fig. |l4|. We observed 
that the resulting changes in the model terminal velocity 
curves are caused in about equal parts by the change in the 
viewing direction relative to the bar, and by the intrinsic 
differences between the deprojected mass distributions. The 
models' scaling factors £ are again determined separately 
for each model by eyeball fitting the observed northern and 
southern terminal curves. The variation in £ is very small, 
however: we obtain £ = 1.1, 1.11, 1.12, 1.11, and 1.12 for 
orientation angles from 10° to 30°. 
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We find only a weak preference for the 20° model. 
The —30° arm is more consistent with the 10° and 15° 
case, whereas these models show too large terminal veloci- 
ties around +40°. The +30° arms seem to be in favour of 
15 — 25° , whereas the 3-kpc-arm, although always too slow 
at I = 0, seems to fit slightly better with 20°. We conclude 
that the orientation is about 20° with a large uncertainty. 



4.6 Spiral arm tangents 



In the previous §4.5 we have already used the observed spiral 



arm tangents at I = ±30° to constrain the pattern speed 
and orientation angle of the bar. Here we reconsider the 
location of the spiral arms in the models that compared best 
with the Galaxy in §4.5 (i.e., those with R c = 3.4 kpc and 
<y2bar = 20°), but with (i) a possible dark halo included, (ii) 
emphasis also on the nearby terminal curve and the spiral 
arms with tangents at I — ±50°, and (iii) higher resolution. 

For comparing the model spiral arms with observations 
we again use the tracers in Fig. |^ (HII regions and molecu- 
lar clouds) and the characteristic features in (l,v) diagrams 
like Fig. oL Unfortunately, because the structures in the ob- 
served (I, v) diagrams are much less sharp than in the model 
(l,v) diagrams, it is very difficult to measure reliably any 
features beyond those already discussed, i.e., the spiral arm 
tangent point positions, the 3-kpc arm and the molecular 
ring. Already the run of the arms out of the molecular ring 
to their tangent points cannot be identified unambiguously 
from Fig. [j]. 

Spiral arm tangent directions in the models are easily 
determined. When the arm is broadened in the tangential 
direction, we place the tangent at the outer edge, where the 
velocity jump is. In this way we can achieve a fair accuracy 
of a few degrees. Only one tangential direction, the Scutum 
arm at I = 30° , cannot easily be determined in this way, be- 
cause its tangent goes through the corotation region through 
which no arm continues inwards in the model. Nevertheless, 
we can measure an approximate value for this tangent as 
well. Table |l| gives a comparison of model and observed spi- 
ral arm tangent directions, and the various data that have 
been used for the observed directions. 

The bar-driven spiral arms in the models change their 
relative strength, opening angle, and hence location, depend- 
ing on pattern speed, bar orient atio n, halo mass contribu- 
tion and other parameters. 



In 



4.5 we have discussed the 



constraints on the the pattern speed and bar orientation an- 
gle. Fig. [li] shows that adding a dark halo mass contribution 
so as to make the outer rotation curve flat (at 208 km s^ 1 
after scaling) has the following two effects: (i) the outer pair 
of spiral arms at I = ±50° now form a more regular pat- 
tern with the inner pair at I = ±30°, and all four arms now 
have a similar density contrast with respect to the interarm 
gas; (ii) all four tangent directions now fit the observations 
reasonably well; see also Table 1. Fig. [iB] also shows the po- 
sitions of large HH-regions and molecular clouds superposed 
on the gas arms. Note that their distances were determined 
from a circular orbit model and so may be slightly in error. 
Nonetheless it is reassuring that they fall approximately on 
the model gas arms, consistent with the model's match to 
the spiral arm tangents. 



4.7 Gravitating spiral arms 

In this section we estimate the effect, if any, of the gravi- 
tational potential of the stellar spiral arms that are likely 
associated with the spiral arms seen in the gas. Remember 
that in our model, the spiral arms outside the bar's corota- 
tion radius are driven by the clumps of NIR light and mass 
3 kpc down the minor axis, which rotate with the bar (see 



§4.1). We interpret these clumps as the signature of real spi- 



ral arms in the deprojected NIR light. In the gas model two 
of the spiral arm heads are about at the correct positions 
where these clumps are observed; this supports the view that 
the gravitational potential of the clumps is a first estimate 
of the true spiral arm potential. 

However, the models discussed sofar did not take into 
account the gravitational perturbations outside the clump 
regions that could be associated with the spiral arms. In 
particular, we are interested in knowing whether the mor- 
phology of the spiral arms would be changed when these 
arms carry a reasonable fraction of the mass throughout the 
disk. To test this, we have experimented with the following 
scheme that makes use of the particular strengths of the SPH 
method: First, we assign a fraction of the mass of the NIR 
disk to the spiral arms. Most of this mass will be associated 
with the spiral arm perturbation in the stellar density. Then 
we assign this part of the mass to the gas particles, assuming 
that a large fraction of the gas will later be concentrated in 
the arms and that the shock fronts seen in the gas will trace 
the stellar spiral arm crests. Finally, to mimic the fact that 
the stellar spiral arms are much broader in azimuth than 
the shocks seen in the gas, we set the gravitational soften- 
ing radius of the gas particles to a value appropriate for the 
stellar arm widths and different from the smoothing length 
used in the calculation of the hydrodynamic forces. 

According to Rix & Zaritsky (1995), the spiral arms 
cover about one third of the surface area in galaxies mor- 
phologically similar to the MW. The arm-interarm contrast 
is somewhat less than average for these galaxies, so that we 
estimate that about 20% of the total NIR luminosity is in the 
spiral arms (superimposed on the axisymmetric background 
disk which thus contributes about 80%). 

The azimuthally averaged surface density associated 
with the spiral arms should thus also be ~ 20% of the mean 
background stellar surface density. In the solar neighbor- 
hood, the stellar surface density of our normalized models is 
35£ 2 M©pc" 2 . Thus we first take an 'arm' potential corre- 
sponding to a mean surface density of 7Mq pc~ 2 at the solar 
radius, and a factor exp[— (r — 8kpc)/2.5 kpc] more further 
in. Here we have used the radial scale-length of the NIR disk 
from §3. This spiral arm mass is given to the gas particles, 
and in order to compensate for the extra mass, we substract 
the same amount as a constant fraction from the COBE NIR 
mass model. 

Because the stellar spiral arms are less sharp than 
gaseaus arms, we estimate their gravitational force by 
smoothing the gravitational forces of the gas particles over 
some length scale e. Each particle contributes a smoothed 
potential 



Mr) 



Gra^i 



(5) 



y/{r — Vi) 2 + e 2 
to the gravitational field of the 'arms'. Here, m*; denotes 
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Figure 15. Spiral arm tangents in two models with y?bar = 20° compared to observations. Left: standard model without dark halo; 
right: including dark halo. Both models are computed with 100,000 particles and the gas disk is truncated at lOkpc. 
The observed directions of spiral arm tangents are shown by the straight lines, starting at the position of the Sun at x = — 7.5kpc, 
y = — 2.7kpc in this plot. See Table |l] which also lists the spiral arm tangents for the models. For illustration, stars denote positions of 
major HII-Regions from Georgelin & Georgelin (1976) and Georgelin et al. (1996), and circles the positions of large molecular clouds from 
Dame et al. (1986) and Grabelsky et al. (1988). The distances of these tracers have been rescaled to Ro = 8kpc but are not corrected 
for non-circular motions. Such a correction would tend to move HII-Regions into the model spiral arms, because of velocity crowding. 



the stellar mass associated with particle i; m,i is obtained 
by dividing the total mass in the stellar spiral arms by the 
effective number of gas particles. The parameter e must de- 
pend on the average distance between two arms in the model 
(about 3kpc in the region of interest). We thus mimic the 
broader spiral arm potential by smoothing the gas arm po- 
tential over about e ~ 1 kpc. 

Notice that the hydrodynamical forces do not depend 
on the actual value of the surface mass density, but only on 
the particle density gradient. We can therefore save some 
memory space by taking the SPH particle mass m; equal to 
the stellar m,i. 

The initial gas disk is set up in a similar way as for the 
models discussed above, but for consistency the surface den- 
sity is that of an exponential disk with the same radial scale 
length as the stellar disk. The resulting additional radial 
pressure gradient is far to weak to change the dynamics. 

In our first attempt at such a model we found that the 



mass in the gas particles accumulating on the S2-disk due to 
inflow was so large that the rotation curve in this region was 
changed significantly. Moreover, the dust lane shock fronts 
acquired significant mass. This is unrealistic since no strong 
stellar arms form at the edge of the bar, and it has the effect 
of changing the gas flow near the cusped orbit. To avoid both 
effects we have set the gas particle masses to zero inside 
the bar region, which was approximated by an ellipse with 
major axis 3 kpc and minor axis 1.8 kpc. Correspondingly, 
the COBE mass model was then not changed in this region. 

The gravitating arm model with the specified parame- 
ters that finally results is shown in the left panels of Fig. |l6| 
First, we notice that the distribution of gas has not changed 
much compared to Fig. |l5| Especially the locations of the 
arms are unchanged. The inner arms now contain more par- 
ticles, because of the initial exponential surface density pro- 
file. However, the line-of-sight velocities have been modified 
somewhat. For example, the 3-kpc-arm has now about the 
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right velocity: ~ — 41kms _1 at I = 0. At around I = 20° 
the terminal velocities are now larger than the observed ve- 
locities. On the southern side, at / ~ —20°, the velocities 
are also higher, removing part of the previous discrepancy 
between model and observation. The larger velocities in this 
region are probably due to the mass in the 3-kpc-arm and 
the second arm at ~ 2 kpc and their symmetric counterarms. 

We have also run a model with about twice as much 
mass in the 'arm' potential and half the smoothing length, 
e = 0.5 kpc (right panels of Fig. |l^). Both parameter in- 
crease the gravitational response to the arms and are ex- 
treme values for the Milky Way. In this case, the distribu- 
tion of the gas has changed more dramatically, especially in 
the region close the the bar. The inner arms seem to form 
an ellipse around the bar; however, the number of arms has 
still not changed. The 3-kpc-arm in this model is faster than 
the observed 3-kpc-arm and expands with 60 kms -1 towards 
the observer at I — 0, but it no longer extends out to 3 kpc. 
The model terminal curve now gives only a poor fit to the 
observations. One interesting property of this model is that 
it contains much more gas with forbidden velocities than all 
our other models. 

From both models it is clear that the gravitational po- 



tential of the spiral arms is important for the comparison 
with Galactic (l,v) diagrams at a level of ~ 15kms . Es- 
pecially the 3-kpc-arm and details in the terminal velocity 
curve depend on this parameter. Besides that, however, the 
morphology and spiral arm tangents are not effected much. 
It is encouraging that the inclusion of spiral arm gravity ap- 
pears to improve the fit to some aspects of the data, even 
though the models discussed in this section certainly do 
not contain the entire story. For example, we have not sub- 
tracted a fraction of the mass in the minor axis NIR clumps 
which should now be taken care of by the spiral arms. On 
a more fundamental level, it is quite possible that the spi- 
ral arms between the bar and the solar circle rotate with a 
somewhat different pattern speed from that of the bar (see 
Sellwood & Wilkinson 1993 for a review of this subject, and 
Amaral & Lepine 1997 for a separate model of the Galac- 
tic spiral arms with a slow pattern speed). If this were the 
case, we would have to observe the resulting time-dependent 
pattern of bar and spiral arms at a moment when both the 
apparent NIR distribution of light and the induced kine- 
matical perturbations resemble those observed in the Milky 
Way. This may not be as difficult as it appears because of 
the Galaxy's relatively tightly wound four-armed spiral pat- 
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tern. However, simulating this would introduce an entirely 
new degree of freedom and will not be attempted here. 



5 CONCLUSIONS 

We have presented a new set of hydrodynamical models 
for the gas flow in the Galactic disk inside the solar circle. 
These gas flows are evolved in the gravitational potentials 
obtained by deprojecting the NIR luminosity distribution of 
the Galactic bar and disk from COBE/DIRBE under the 
assumption of eight-fold symmetry (as in Binney, Gerhard 

6 Spergel 1997), assuming constant mass-to-light ratio for 
the NIR luminous material, and adding a nucleus and (in 
some cases) a dark halo component. These models allow us 
to understand many features of Galactic HI and CO obser- 
vations. 

To follow the gas dynamics we have used the Smooth 
Particle Hydrodynamics (SPH) code described in Englmaier 
& Gerhard (1997). With the SPH method we can resolve 
the spiral arm shocks well, and use them in some models 
as tracers for the gravitational potential of the stellar spi- 
ral arms. The hydrodynamical models quickly settle to an 
approximate quasi-equilibrium flow pattern whose overall 
morphology is not sensitive to the precise value of the bar's 
pattern speed, to the orientation of the bar with respect to 
the observer, and to whether or not the spiral arms carry 
mass. 

We have compared our gas models with Galactic HI, 
CO, HII, and other data. We find that these models provide 
a coherent explanation of many aspects of the data, such 
as: (i) the four-armed spiral structure of the Milky Way 
between corotation and the solar radius, (ii) the nature of 
the 3-kpc-arm, (iii) the terminal velocity curve, (iv) the 
non-circular velocities near the cusped orbit at the ILR, 
and (v) the disk of gas on the inner X2-orbits. Thus NIR 
photometry and gas kinematic observations conform to a 
single picture, and the Galactic bar is an essential part of 
this. 

In this picture, the bar (bulge) rotates with a pattern 
speed such that corotation is at R c ~ 3.5 ± 0.5 kpc. The 
3-kpc arm is one of the arms emanating from the ends of 
the bar, extending into the corotation region. Outside coro- 
tation, a four-armed spiral arm pattern gives rise to the 
molecular ring and the arms extending to the solar circle 
and beyond. In the model, this pattern is generated by the 
rotating luminosity/mass concentrations on the bar's minor 
axis found by BGS. These can therefore not just be due 
to light from young supergiant stars, but must be massive; 
most likely they are symmetrized approximations to the stel- 
lar spiral arms themselves. A spiral pattern similar to that 
found here has been observed by Fux (1998) in those of his 
N-body - SPH barred galaxy models which compare best to 
the Galactic (l,v) diagrams. These models start from a set 
of specified initial conditions rather than from observations. 
The fact that both approaches lead to similar overall results 
is encouraging. 

We find that the Galactic terminal curve out to longi- 
tudes I ~ 45° is consistent with a maximal, constant mass- 
to-NIR light disk and bar model. The inferred mass in the 
disk cannot easily be reduced because (i) this model still un- 
derpredicts the microlensing optical depth towards the bulge 



(Bissantz et al. 1997) and (ii) it predicts about the correct 
surface mass density for the old stellar disk near the Sun. 
Thus the Galactic dark halo will be an important contribu- 
tion to the mass of the Milky Way only outside of at least 
R — 5 kpc, depending on the LSR rotation velocity. 

The models are similar and in reasonable overall agree- 
ment with the Milky Way observations for a range of val- 
ues of the bar's orientation angle. Probably best is t^bar = 
20 — 25° but the uncertainties are are such that 15° or 30° 
cannot be excluded. The detailed match to various key ob- 
servational quantities does depend on ipb&r and the precise 
value of R c , but at a level comparable to the influence of 
various other parameters including the LSR rotation veloc- 
ity, the asymptotic halo circular speed, and the influence of 
the stellar spiral arm gravity. 

Some features of the HI and molecular emission data 
can be reproduced quantitatively, such as the Galactic spi- 
ral arm tangent positions, much of the terminal velocity 
curve, or the position of the nuclear ^2-disk. On the other 
hand, a detailed quantitative fit to all features in the ob- 
served (I, v) diagrams is not yet possible. The 3-kpc arm 
in the model has the correct angular extent but somewhat 
too small non-circular velocities. Its arm tangent position 
coincidences with the place where the model terminal curve 
differs most from the observations. In this region the disk 
mass model may not be very accurate. The position and ve- 
locity of the molecular parallelogram are fitted well by closed 
orbits in the NIR potential, whereas the hydrodynamic gas 
flow underestimates the velocity and overestimates its radial 
scale. Probably both hydrodynamic resolution and uncer- 
tainties in the potential near ~ 200 pc (where the NIR data 
have insufficient resolution) are responsible for this. Also, 
this suggests that the clouds near the peak of the terminal 
velocity curve have a small mean free path. 

Much work remains to be done. Because the deprojected 
disk model accounts only insufficiently for the Galactic spiral 
arms, some aspects of the gravitational potentials used are 
likely to be wrong. With a model for the spiral arms in 
hand, both the correction of the NIR data for dust and the 
subsequent deprojection of these data could be improved. 
Further observational work on spiral arm tracers such as 
HII regions and molecular clouds would be highly valuable 
for clarifying the run of the weaker Galactic spiral arms 
and those on the other side of the bar, and thus for better 
constraining the gas-dynamical models. Distance estimates 
to these tracers could be improved by making use of the 
velocity fields in these models. 

On a more fundamental level, several assumptions made 
in the models may be or are likely to be invalid at some 
level, and require further study: (i) that the gravitational 
potential and gas flow structures are quasi-stationary and 
point-symmetric with respect to the Galactic Centre; (ii) 
that the NIR light is a fair tracer of the stellar mass, i.e., 
that young supergiant stars do not contribute significantly 
to the NIR light (but see Rhoads 1998); (iii) that the gas 
disk can be treated as planar; (iv) and that the bar and 
spiral arms rotate with the same pattern speed. 
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6 MULTIPOLE EXPANSION 

For the SPH simulations, the gravitational forces are cal- 
culated using a multipole expansion of the stellar density 
model. The stellar density is given by 



p(r,6,ip) = pof(r,0,<p) 



(6) 



where the function / is the interpolated deprojected light 
distribution within the bulge region and the exponential 
disk elsewhere. Here, the constant po is the presently un- 
known unit of luminosity density multiplied by the also 
unknown mass-to-light ratio. A first estimate gives po = 
£ 2 3 10 8 Mq/pc 3 , where £ is of the order of unity and is 
measured by fitting the observed rotation velocity in § 4.4 . 
The density multipoles 

/• TT /•27T 

Pi m (r) = / d8 sm(9)P lnl (cos(9)) / dip cos(rmp) f (r) (7) 



satisfy the identity equation 

(2l + l)(l-m)\ 

I . rri 



p(r) = ^(2-<5 m0 ) 



47r(7 + ' 



PlmPlmCOsijrup). (8) 



The Pi m are the associated Legendre functions. Here we have 
used that / is an even function in <j>: f(r, 6, 4>) = f(r, 0, ~<f>), 
and have restricted the sum to positive m. 

Both integrations were performed with the Romberg 
method. The density multipoles pio were then fitted to a 
power law Cir Pl in the range of 350 pc to 500 pc with the 
method of least squares. The tabulated multipole expansions 
were then replaced by the fit inside 350 pc. 

From the modified tables pi m (r) we calculated the fol- 
lowing two auxiliary integrals: 



I<( r ) = / dap tm (a)a 
Jo 

and 



1 + 2 



J>(r) = 



da pi 

in (a) a 



(9a) 



(9b) 



For stability reasons we used for these integrations the trape- 
zoidal rule on the 3000 logarithmically equidistant pi m val- 
ues tabulated between 1 pc and 12 kpc. The region between 
12 and 16 kpc in the second integral was calculated again 
with the Romberg method; outside 16 kpc the density was 
set to zero. 

From this, we get the potential multipoles 

$ im (r) = (2-5 m0 )|^|(r- i - 1 /<(r)+r ! 7 > (r)) (10) 



as well as their first derivatives 

*' lra (r) = (2-o m0 )| T ^||(-(/+l)r- i - 2 / < (r)+Z r'" 1 /> (r)).(ll) 
The potential is then given by 

$(r) = -Gp ^2 ®im(r)Pi m cos(rrvp). (12) 

l,m 

The components of the gravitational acceleration can 
be calculated from the partial derivatives of $, 

d$ 



— = -Gp ^2 ^ImPlm COs(mip), 
l,m 

^| = Gpo— ^ ®lmPlm COs(nvp), 

do r 

l.m 

^— = Gpo m ®imPim sin(m</p), 



and the components of the acceleration a = — X7& from 



d$ d$ 2 d$ 2 

a v = --rvn - i- x / R - -TEV z l\ r R )' 



(13) 
(14) 
(15) 

i 

(16) 
(17) 

dr~" ' 6B~>- ' (18) 
For the gasdynamical model we need only the forces in 
the galactic plane. In this special case, it is therefore suffi- 
cient to tabulate the functions 



dr ' dip 



de- 



d$ d$ 2 

a z = -—z/r + -T^R/r 



$m = 2J ®lmPlm(0)$lm, 
I 

*m = ^2&lmP m (,Q)<$>l m , 
I 

and to compute the forces using 
— = -Gpo ^2 cos(m^), 

m 

^— = Gpo m$ m sinOrrup). 
dip * — ' 



(19) 
(20) 

(21) 
(22) 
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